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SUMMARY 


"Limit  Design"  guides  the  structural  engineer  towards  an 
economic  design  of  structures  that  are  made  of  ductile  materials 
and  have  to  carry  specified  loads.  Prom  the  mathematical  point 
of  view,  the  problem  can  be  reduced  to  one  in  linear  programming, 
but  even  simple  structures  may  lead  to  linear  prografnming  prob¬ 
lems  of  considerable  size.  A  method  of  solution  is  discussed 
that  has  been  found  efficient  for  structures  of  moderate  com¬ 
plexity.  A  sketch  of  the  historical  development  of  limit  analy— 

» 

sis  and  design  is  given. 
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LINEAR  PROGRAMMING  AND  STRUCTURAL  DESIGN 
II .  Limit  Design* 

William  Prager 

In  general  terms,  the  problem  discussed  in  the  preceding 

** 

lecture  can  be  described  as  follows;  given  the  geometry  of  a 
structure  and  the  fully  plastic  moments  of  its  members,  to  find 
its  carrying  capacity  for  a  given  type  of  loading.  This  is 
known  as  the  problem  of  "limit  analysis."  Today's  lecture  will 
be  primarily  concerned  with  the  related  but  different  problem 
of  "limit  design":  given  the  geometry  of  a  structure  and  the 
loads  it  has  to  carry,  to  choose  the  fully  plastic  moments  of 
its  members  so  as  to  minimize  the  weight  of  the  stiMcture. 

Before  entering  on  the  discussion  of  limit  design,  let  us 
review  the  method  developed  last  time.  Figure  la  shows  a  portal 
frame;  the  feet  of  the  columns  are  build  in,  and  the  beam  is 
rigidly  Joined  to  the  colxomns.  The  fully  plastic  moment  of 
each  column  is  Mq  and  that  of  the  beam  is  2Mq.  For  the  type 
of  loading  Indicated  in  the  figure,  we  wish  to  determine  the 
maximum  load  Intensity  that  the  frame  can  carry;. 

I 

i 

The  potentially  dangerous  cross  sections  of.  the  frame  are 
at  the  foot  of  each  column,  at  the  top  of  each  column  Just  below 
the  Joint  with  the  beam,  and  at  the  points  of  application  of  the 
loads.  These  cross  sections  are  numbered  1  through  6  In  Fig.  la, 

*Iiecture  delivered  at  RAND  on  July  11,  1957. 

Linear  Programming  and  Structural  Design,  I.  Limit  Analysis 
The  HARD  Gorporation,  Paper  P-1122,  June  2b,  1957 • 
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and  the  bending  moments  at  these  sections  will  be  denoted  by 
through  These  bending  moments  must  satisfy  three  equa¬ 

tions  of  equilibrium,  which  may  be  established  kinematically 
by  considering  the  three  "basic  mechanisms"  Indicated  In  Pigs, 
lb,  Ic,  and  Id.  The  coefficients  and  right-hand  sides  of  these 
equations  are  listed  In  rows  (a)  through  (c)  of  Table  I,  bend¬ 
ing  moments  being  considered  positive  when  they  tend  to  produce 
tensile  stresses  in  the  fibers  bordering  on  the  Interior  of  the 
frame . 

TABLE  I 


t  ”2 

"j 

«4 

"5 

«6 

r.h.s . 

Pa/ho 

(b) 

5 

-2 

0 

B 

IB 
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5.00 

(c) 

0 

0 

B 

2 

B 

B 

2Pa 

5.00 

(d) 

-1 

0 

1 

0 

B 

B 

2Pa 

2.00 

«i/«0 

-1 

2,67 

B 

2 

B 

1 

— 
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(e) 

0 

5 

-5 

0 

B 

B 

0 

— 

(f) 

-5 

5 

B 

0 

-2 

2 

oPa 

1.67 

«iAo 

-1 

1 

•5^ 

1 

1.55 

-1 

1 

— 

— 
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Table  I  contains  all  the  information  necessary  for  for¬ 
mulating  the  linear  programming  problem  that  furnishes  the 
load  carrying  capacity  of  the  frame.  We  may,  for  instance, 
take  the  Equation  (b)  as  defining  our  objective  function  and 
obtain  two  subsidiary  conditions  for  the  unknown  bending  mo¬ 
ments  by  eliminating  P  from  Equations  (b)  and  (d).  We  thus 

t 

obtain  the  following  problem;  to  maximize 

2Pa  -  +  3M2  -  2M^,  (1)  . 

subject  to  the  equations 

Ml  -  3M2  +  +  2M^  -  -  0,  (2) 

3M2  +  3M5  -  M5  +  Mg  -  0,  (3) 

and  the  inequalities  stipulating  that  Mj^  must  not  exceed  2Mq 
in  absolute  value  while  the  other  five  bending  moments  must 
not  exceed  Mq  in  absolute  value. 

If  we  have  a  computer  code  for  linear  programming  with 
bounded  variables,  we  can  at  this  stage  put  the  problem  on 
the  computer  and  have  it  compute  the  load  carrying  capacity  of 
the  frame.  In  the  present  case,  this  would  hardly  be  worth 

t 

while  because  the  problem  is  readily  solved  by  the  method 
developed  in  the  last  lecture.  Considering  in  turn  each  of 
the  mechanisms  of  Pigs,  lb  through  Id  as  a  tentative  collapse 
mechanism,  we  give  the  appropriate  values  to  the  bending  mo¬ 
ments  at  the  plastic  hinges  and  compute  Pa/M^  by  applying  the 
principle  of  virtual  work;  the  results  are  indicated  in  the 
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last  column  of  Table  I.  The  mechanism  of  Pig.  Id  is  found  to 
give  the  lowest  value,  Pa/Mg  =  2,  the  bending  moments  at  the 
plastic  hinges  being  »  -Wq  and  «  Mg  *  Mq.  Inserting 

these  values  into . Equations  (b)  and  (c),  we  solve  for  Mg  and 
M|^,  the  bending  momenta  at  the  sections  where  the  considered 
mechanism  does  not  have  plastic  hinges.  The  row  just  below 
(d)  in  Table  I  lists  the  resulting  dimensionless  bending  mo¬ 
ments  M^/Mq.  As  Mg  exceeds  its  upper  bound  Mq,  we  must  super¬ 
impose,  on  the  considered  mechanism  (Pig.  Id),  another  mech¬ 
anism  that  has  a  positive  hinge  at  2  and  no  hinge  at  4,  and 
yields  zero  total  work  of  the  loads.  This  mechanism  is,  of 
course,  a  linear  combination  of  the  three  basic  mechanisms 

V 

(Pigs,  lb  through  Id).  To  within  an  Irrelevant  factor,  the 
coefficients  in  this  linear  combination  are  determined  by  the 
following  conditions;  in  the  equation  of  equilibrium  cor¬ 
responding  to  the  desired  mechanism,  the  coefficient  of  Mg 
must  be  positive,  while  that  of  and  the  right-hand  side 
must  vanish.  A  glance  at  Table  I  shows  that  the  combinaition 
(b)  —  (d)  satisfies  these  conditions.  The  resulting  equation 
of  equilibrium  is  given  in  row  (e),  and  the  corresponding 
mechaunism  is  shown  J.n  Pig,_le.  Comparing  the  coefficients 
in  the  rows  (d)  and  (e),  we  note  that  the  smallest  positive 
value  A  for  which  a  coefficient  in  the  combination  (d)  +  X(e) 
vanishes  is  A  ■  l/5.  To  avoid  fractions,  we  consider  the 
combination  5(d)  +  (e)  which  is  given  in  row  (f).  The  cor-  '  ■ 
responding  value  of  Ps/Mq  is  I.67  as  indicated  in  the  last 
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column  of  row  (f).  When  the  bending  moments  are  computed  as 
before,  the  values  in  the  last  row  of  Table  I  are  obtained. 

Since  these  are  all  within  the  bounds  for  the  respective  mo¬ 
ments,  we  have  found  the  load  carrying  capacity  of  the  frame 
to  be  P  »  1,67  Mo/a.  The  corresponding  collapse  mechanism  is 
shown  in  Pig.  If. 

At  this  stage  it  seems  appropriate  to  discuss  several 
questions  which,  I  am  sure,  have  occurred  to  some  of  you. 

First  of  all,  is  this  systematic  approach  really  necessary, 
or  would  the  collapse  mechanism  be  obvious  to  an  experienced 
analyst?  For  the  simple  problem  considered  here,  one  would 
certainly  consider  the  mechanism  of  Pig.  If  as  a  likely  collapse 

4 

mechanism  because  it  does  not  have  plastic  hinges  in  the  rela¬ 
tively  strong  beam.  Computing  first  the  load  intensity  neces¬ 
sary  to  activate  this  mechsinism,  and  then  the  bending  moments 
at  the  sections  2  and  4  as  before,  one  easily  verifies  that 
this  is  indeed  a  collapse  mechanism  of  the  frame.  In  more 
complex  structures,  however,  the  initial  guess  regarding  the 
collapse  mechanism  is  likely  to  prove  erroneous;  the  systema¬ 
tic  procedure  described  above  will  then  lead  to  the  correct 
mechanism.  i 

I 

A  question  that  forces  Itself  on  anybody  familiar  with 
linear  programming  concerns  "degeneracy."  Described  in  struc¬ 
tural  terms,  degeneracy  occurs  when  the  mechanism  obtained  at 
a  certain  stage  of  our  procedure  has  less  than  the  "regular", 

I 

number  of  plastic  hinges.  For  the  frame  of  Pig.  la,  for 
instance,  we  have  six  possible  locations  of  plastic  hinges 
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and,  after  eliminating  P  between  the  three  equilibrium  con¬ 
ditions,  two  linear  homogeneous  equations  for  the  bending  mo¬ 
ments  through  M^.  These  allow  us  to  compute  two  bending 
moments,  when  the  four  others  have  been  chosen  in  accordance 
with  the  relative  rotations  in  the  four  plastic  hinges  of  a 
tentative  collapse  mechanism.  For  the  frame  of  Pig.  la,  the 

f 

mechanisms  in  Pigs.  Id  through  If,  which  have  four  hinges, 
are  therefore  "regular,"  while  those  in  Pigs,  lb  and  Ic,  which 
have  only  three  hinges,  are  "degenerate." 

In  the  theory  of  linear  programming,  degeneracy  is  es¬ 
sentially  treated  by  slightly  modifying  the  problem  so  as  to 
avoid  degeneracy  in  the  modified  problem.  To  show  that  a 
similar  treatiment  of  degeneracy  may  be  possible  in  limit  anal¬ 
ysis  of  structures,  let  us  again  consider  the  frame  in  Pig,  la, 
assuming  that  the  tentative  collapse  mechanism  is  chosen  ac¬ 
cording  to  Pig.  lb.  The  three  hinges  of  this  mechanism  allow 
Infinitesimal  horizontal  displacements  of  the  cross  sections 
of  the  left-4iand  column  only  because  they  lie  on  the  same 
straight  line.  It  is  therefore  natural  to  consider  the  frame 
of  Fig.  la  as  the  limiting  case  of  frames  that  have  slight 
built-in  kinks  at  2  (Pig.  ^).  To  avoid  fractional  values  of 

j 

linear  or  angular  displacements,  the  depth  of  the  kink  is 
denoted  by  4€a  in  Pig.  2.  In  Fig.  lb,  the  segments  1,  2  and 
2,  5  of  the  left-hand  column  have  rotated  by  the  small  amounts 
9  and  20  in  the  clockwise  and  counterclockwise  senses,  re¬ 
spectively.  For  the  frame  of  Pig.  2,  the'  same  rotations  would 
produce  a  vertical  downward  displacement  of  point  5  equal  to 


» 


P-1125 

7-11-57 

-7- 


12ea0,  the  horizontal  displacement  of  this  point  being  zero 
as  in  Pig.  lb.  To  make  this  vertical  displacement  possible, 
we  must  insert  an  additional  plastic  hinge  at  either  4  or  5 . 

The  first  choice  would  not  alter  the  load  intensity  P  ■  5MQ/a 
needed  to  activate  the  mechanism,  whereas  the  second  choice 
lowers  it  to  5M^[(l+3e)a]  .  We  therefore  choose  the  additional 
hinge  at  5,  setting  =  -Mq  and  Mg  «  Mq.  Proceeding 

as  before,  we  are  immediately  led  to  the  collapse  mechanism 
of  Fig.  If.  While  it  seems  likely  that  degeneracy  can  always 
be  treated  in  this  fashion,  no  general  proof  of  this  seems  to 
have  been  given  so  far. 

Finally,  there  is  the  question  of  whether  it  is  realistic 
to  restrict  the  discussion  to  a  single  "type  of  loading."  With 
reference  to  Pig.  la,  for  Instance,  we  might  assume  the  vertical 
and  horizontal  loads  to  vary  independently  between  given  limits, 
and  ask  whether  any  of  the  resulting  states  of  loading  over— 
teixes  the  frame.  Denoting  the  horizontal  load  by  P  as  before, 
and  the  vertical  load  by  Q,  we  then  have  to  consider  all  pos¬ 
sible  collapse  mechanisms  and  specify  that  for  each  of  them 
the  work  dissipated  in  the  plastic  hinges  should  not  be  less 
tHan  the  work  of  the  applied  loads.  For  the  mechanism  of 
Pig.  Id,  for  Instance,  this  condition  provides  the  Inequality 

2Mq  >  |Pa|  ,  (4) 

where  the  absolute  value  has  been  used  on  the  right-hand  side 
because  negative  values  of  P,  l.e.,  horizontal  loads  directed 
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towards  the  left  In  Pig.  1,  may  be  admitted.  Using  Ps/Mq  and 
Oa/MQ  as  rectangular  Cartesian  coordinates  of  the  "load  point" 
in  a  "load  plane,"  we  see  that  (4)  restricts  the  load  point 
to  a  strip  that  is  centered  at  the  origin.  Each  potential 
collapse  mechanism  furnishes  such  a  strip;  the  polygon  common 
to  all  strips  (Pig.  3)  contains  the  load  points  corresponding 

f 

to  states  of  loading  within  the  carrying  capacity  of  the  frame . 
On  the  other  hand,  the  given  ranges  of  variation  of  P  and  Q 
specify  a  rectangle  in  the  load  plane  and  are  admissible  only 
if  this  rectangle  is  contained  within  the  polygon  of  Pig.  3. 

Let  us  now  turn  to  the  problem  of  limit  design:  given 
the  geometry  of  a  structure  eund  the  loads  it  has  to  carry, 
to  determine 'the  fully  plastic  moments  of  its  members  so  as 
to  minimize  its  weight.  Por  the  treatment  of  this  problem 
it  is  customary  to  assume  that  the  weight  per  unit  length  of 
a  structural  member  is  proportional  to  its  fully  plastic  mo¬ 
ment  Mq.  This  assumption  greatly  simplifies  the  einalysis 
because  it  leads  to  a  linear  rather  than  a  nonlinear  programming 
problem.  Whether  the  assumption  is  reasonable  in  a  particular 
case,  depends  on  the  range  of  cross-sectional  shapes  that  are 
at  the  disposal  of  the  designer.  If,  for  Instance,  all  avail¬ 
able  cross-sections  are  geometrically  similar,  the  unit  weight 

a/3 

is  proportional  to  Nq  ,  as  a  simple  dimensional  analysis  will  < 
show.  Actually,  the  standard  I  sections  are  not  strictly 
similar  to  each  other,  and  the  exponent  2/5  must  be  replaced 
by  one  closer  to  unity.  Por  a  more  detailed  discussion  of 


t 
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this  point,  Heyman's  [l]*  book  may  be  oonsulted.  In  the  fol¬ 
lowing,  the  weight  per  unit  length  of  a  member  will  be  taken 
as  proportional  to  its  fully  plastic  moment. 

Figure  4a  shows  a  beam  on  three  supports  subjected  to 
specified  vertical  loads.  The  unknown  fully  plastic  moments 

of  the  left—  and  rlght^and  spans  will  be  denoted  by  Mq  and 

*  ’’ 

Mq,  respectively.  Any  plastic  hinge  at  the  central  support 

will  form  Just  to  the  left  or  right  of  this  support  depending 

*  *  • 

on  whether  Mq  <  Mq  or  Mq  <  Mq.  Since  it  is  not  known  before¬ 
hand  which  of  these  inequalities  api^lies,  we  must  consider 

% 

the  four  potentially  dangerous  cross  sections  marked  1  to  4  in 
Pig.  4a.  The  equilibrium  conditions  for  the  bending  moments 
M^  through  M||  at  these  sections  may  be  derived  from  the  mechan¬ 
isms  shown  in  Pigs.  4b  through  4d.  The  first  two  of  these 
figures  show  the  familiar  "beam  mechanisms,"  the  last  a  "Joint 
mechanism;"  the  infinitesimal  length  of  beam  between  the 
sections  2  emd  3  is  tilted  to  produce  finite  rotations  in  the 
hinges  2  and  3  but  only  infinitesimal  deflections.  The  equi¬ 
librium  conditions  for  the  four  bending  moments  are 


2M^  -  Mg  -  50, 

1 

1 

i 

(5) 

-  +  3Mi^  -  20, 

1 

1 

(6) 

-  M,  -  0 

(7) 

where  the  coefficients  on  the  left-hand  sides  are  dimensionless. 


Numbers  in  square  brackets  refer  to  the  Bibliography  at 
the  end  of  this  report.  -  - 
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whlle  the  right-hand  sides  represent  foot-tons. 

Considering  a  set  of  values  through  Mj^  that  satisfy 
(5) I  (6),  and  (7),  choose 

Mq  -  max  (|Mj|,  IMgl)  (8) 

and 

M*  =  max  (|Mj|,  _  (9) 

The  beam  with  these  fully  plastic  moments  cannot  fall  under 
loads  that  have  the  same  ratio  as  in  Pig.  4a  but  lesser  in¬ 
tensities.  This  beam  therefore  represents  a  feasible  design, 

» 

but  does  not  necessarily  minimize  the  weight,  which  la  pro¬ 
portional  to 

'  W  »  20Mq  +  30M*  .  (10) 

Our  problem  is  to  minimize  W  as  given  by  (10),  where 
Mq  and  Mq  are  defined  by  (8)  and  (9)  in  terms  of  four  bending 
momenta  that  must  satisfy  (5)  through  (7).  This  can  be  for¬ 
mulated  as  a  conventional  problem  in  linear  programming  by 
observing  that  (8)  is  equivalent  to 

Mq  -  Ml  >  0,  j  Mq  ±  Mg  >  0,  (ll) 

! 

} 

and  that  similar  inequalities  niay  be  used  to  replace  (9).  It 
is  worth  noting,  however,  that  even  our  extremely  simple  struc¬ 
tural  problem  leads  to  a  linear  programming  problem  Involving 
6  unknowns  subject  to  thz'ee  equations  and  8  inequalities. 

As  is  to  be  expected,  this  static  formulation  of  the 
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problem  of  mimimum  weight  design  has  Its  kinematic  counterpart . 
The  unknowns  in  this  are  the  fully  plastic  moments  Mq  and  MqJ 
the  conditions  which  they  must  satisfy  are  derived  from  the 
principle  that  for  any  tentative  collapse  mechanism  the  work 
dissipated  in  the  plastic  hinges  should  not  be  less  than  the 
work  done  by  the  loads . 

f  • 

For  the  beam  in  Pig.  4a  all  mechanisms  can  be  considered 
as  linear  combinations  of  the  basic  mechanisms  shown  in  Figs. 

4b  through  4d.  A  tentative  collapse  mechanism  must  have  two 
of  the  four  possible  plastic  hinges. inactive .  Table  II  shows 
the  equations  of  equilibrium  corresponding  to  the  six  potential 
collapse  mechanisms ;  the  first  three  rows  correspond  to  the 
mechanisms  of  Pigs.  4b,  4c,  and  4d,  and  the  last  three  rows 
to  combinations  with  two  inactive  hinges. 

TABLE  II 

II  Mj  Mg  [  ft.t. 

_ (M 

_ _ ^ 

j ^ 

;  (e)  »  (b)  +  (d) 

(f)  -  (c)  -  (d) 

(g)  »  (b)  -  (c)  +  (d) 

fully  plastic  moment 
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The  hinge  rotations  listed  in  coliomns  and  go  with 
the  fully  plastic  moment  Mq  and  those  in  columns  and 
with  Mq  as  indicated  at  the  bottom  of  Table  II.  Since  the 
sign  of  the  bending  moment  at  a  plastic  hinge  corresponds  to 
that  of  the  hinge  rotation,  the  work  dissipated  in  the  hinges 
of,  say,  the  left-hand  span  equals  the  product  of  Mq  by  the 

t 

sum  of  the  absolute  values  of  the  hinge  rotations  of  this  span . 
The  mechanisms  of  Table  II  therefore  yield  the  following  in¬ 
equalities  ; 


3Mq  >  30  , 

(b) 

4Mj  >  20  , 

(c) 

Mq  +  M  J  >  0  , 

(d) 

2Mq  +  M;  2  50  , 

(e) 

Mq  +  3MJ  2  20  ■, 

(f) 

2Mq  +  3MJ  2  10  . 

(g) 

The  problem  of  minimizing  W  as  given  by  (10)  subject  to 
the  inequalities  (b)  through  (g)  is  again  a  problem  in  linear 
programming.  It  is  convenient  to  interpret  these  inequalities 
geometrically  by  using  Mq  and  Mq  as  the  rectangular  Cartesian 
coordinates  of  the  "design  point"  in  the  "design  plane"  (Pig. 
5).  Each  of  the  inequalities  restricts  the  design  point  to  a 
half-plane;  the  "region  of  admissible  designs"  that  is  common 
to  these  half-planes  is  shaded  in  Fig.  3.  According  to  (10), 
all  design  points  that  correspond  to  a  given  structural  weight 
lie  on  a  straight  line,  e.g.,  the  dotted  line  in  Fig.  5*  Since 
this  line  does  not  contain  a  point  of  the  shaded  region,  no 
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deaign  of  this  structural  weight  can  support  the  given  loads . 

As  the  lines  corresponding  to  various  fixed  values  of  W  are 
parallel  to  each  other,  the  desired  design  of  minimiam  weight 
is  obviously  represented  by  the  vertex  V  of  the  region  of 
admissible  designs .  Since  each  of  the  sides  through  this  vertex 
correspond  to  a  collapse  mechanism,  the  minimum  weight  design 
admits  two  competing  collapse  mechanisms.  We  observe  that  the 
number  of  competing  collapse  mechanisms  equals  the  number  of 
dimensions  of  our  "design  space."  When  n  fully  plastic  momenta 
are  at  the  choice  of  the  designer,  the  minimum  weight  design 
must  admit  at  least  n  competing  collapse  mechanisms,  but  may 
have  more  than  n  in  special  cases. 

Making  A  guess  concerning  the  competing  collapse  mechan¬ 
isms  for  the  minimum  weight  design,  one  may  compute  bending 
moments  to  ascertain  whether  this  guess  corresponds  to  a  vertex 
of  the  region  of  admissible  designs.  If,  for  instance,  we 
take  (b)  and  (f)  as  competing  mechanisms  in  our  example,  we 
find  (in  ft.t.)  Mq  ■  10  and  Mq  »  3*35 •  The  bending  moments 
at  the  plastic  hinges  in  either  mechanism  must  then  equal  the 
appropriate  fully  plastic  moments  in  absolute  value  and  have 

t 

the  same  signs  as  the  hinge  rotations.  We  thus  have  (lir  ft.t.) 

Mj  --Mj  ■  10  and  ■  3.33;  equilibrium  finally  requires  that 

* 

Mj  ■  —10,  and  this  exceed  the  fully  plastic  moment  Mq  in 
absolute  value.  Accordingly,  the  assumed  combination  of  col¬ 
lapse  mechanisms  does  not  correspond  to  a  vertex  of  the  regJ^ 
of  admissible  designs.  Applying  the  same  technique  to  the 


« 
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combinatlon  of  (c)  and  (e),  we  find  (In  ft.t.)  Mq  =  12.5  and 

* 

Mq  ■  5.0.  Ihe  bending  momenta  with  plastic  hinges  are  ■  -12.5 
and  M^  =  «  5.0,  and  equilibrium  requires  that  Mg  »  -5.0, 

which  is  less  in  absolute  value  than  Mq.  This  shows  that  we 
now  have  a  vertex  of  the  region  of  admissible  designs,  but 
does  not  reveal  whether  we  have  attained  the  minimum  weight. 

The  computation  of  bending  moments  that  are  compatible 
with  both  collapse  mechanisms  incidentally  yields  the  following 
result.  If  the  competing  mechanisms  have  a  plastic  hinge  in 
common,  their  rotations  at  this  hinge  must  have  the  same  sign. 

Any  possible  linear  combination  of  the  inequalities  (c) 
and  (e)  restricts  the  design  point  to  a  half-plane  that  is 
bounded  by  a  tllne  through  V  and  contains  the  region  of  admis¬ 
sible  designs.  Since  the  line  of  min  W  is  of  this  type,  it 
must  be  possible  to  combine  the  inequalities  (c)  and  (e)  with 

positive  factors  in  such  a  manner  that  the  coefficients  of 

« 

Mq  emd  Mq  in  the  resulting  inequality  coincide  with  those  in 
the  expression  for  W  (Equation  10).  As  has  Just  been  shown, 
the  corresponding  combination  of  the  mechanisms  (c)  and  (e) 
cannot  bring  about  any  partial  or  comp}.ete  cancellation  of 
rotations  in  the  common  hinge  5*  Accordingly,  the  inequality 

i 

obtained  as  a  positive  linear  combination  of  the  inequalities 
(c)  and  (e)  directly  expresses  the  fact  that  in  the  correspond¬ 
ing  combination  of  the  mechanisms  (c)  and  (e)  the  work  dis¬ 
sipated  in  the  plastic  hinges  must  not  be  smaller  than  the  work 
done  by  the  loads.  The  minimum  weight  design  therefore  adnd'ta 
a  collapse  mechanism  for  which  the  sums  of  the  absolute  values 
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of  the  hinge  rotations  associated  with  each  fully  plastic  mo¬ 
ment  have  the  same  ratios  as  the  coefficients  of  these  moments 
in  the  expression  for  the  structural  weight.  In  our  example, 
this  mechanism  is  given  by  5(c)  +  10(e);  it  Involves  the  hinge 

rotations  20  at  1,  —15  at  3,  and  15  at  4. 

» 

For  the  beam  in  Pig.  4  it  was  an  easy  matter  to  list  all 

f 

potential  collapse  mechanisms  and  the  corresponding  inequalities 
For  a  more  complex  structure,  however,  this  may  be  a  formidable 
task.  Once  this  is  accomplished,  the  problem  is  readily  set 
up  as  a  problem  in  linear  programming,  which  may  be  put  on  a 
computer.  Unfortunately,  no  systematic  method  has  as  yet  been 
discovered  for  reducing  the  labor  of  listing  all  potential  col¬ 
lapse  mechani^sms.  The  criteria  developed  above,  however,  are 
useful  aides  in  the  search  for  the  proper  combination  of  mecha¬ 
nisms.  To  Illustrate  this,  consider  the  two-storey  frame  of 
Pig.  6a.  The  feet  of  the  columns  are  built  in,  and  all  Joints 
are  rigid.  There  are  11  potentially  dangerous  cross  sections, 
labelled  1  through  11  in  Pig.  6a.  The  bending  moments  at  these 
sections  will  be  denoted  by  through  Three  fully  plastic 

moments  are  at  the  choice  of  the  designer;  Mq  for  the  columns 
of  the  first  storey,  Mq  for  the  ceiling  beam  of  this  storey, 
and  for  the  columns  cind  the  ceiling  beam  of  the  second 
storey.  Accordingly,  we  must  look  for  three  competing  collapse 
mechanisms.  Five  basic  mechanisms  are  shown — ^Pigs.  6b  through 
6f;  the  corresponding  equations  of  equilibrium  are  listed  in 
Table  III. 


« 
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The  structural  weight  is  proportional  to 

W  -  12Mq  +  8Mq  +  20Mq  (12) 

The  competing  mechanisms  are  likely  to  Include  (b)  and  (c), 
the  first  because  It  determines  the  minimum  fully  plastic  mo¬ 
ment  of  the  upper  storey,  and  the  second  because  It ^requires 
a  particularly  large  amotint  of  work  to  be  dissipated  In  the 
plastic  hinges .  The  Inequalities  corresponding  to  these  mecha—  ' 
nlsms  are 

4Mq2^0  ,  (15) 

and 

j  2Mq  +  4Mq  +  2m”  >  520  .  (14) 

As  (d)  furnishes  the  Inequality 

0 

4M^  >  140,  (15) 

It  would  not  be  possible  to  obtain  the  coefficients  In  (12)  by 
a  positive  linear  combination  of  the  coefficients  In  (15) 
through  (15).  Thus,  If  the  competing  collapse  mechanisms  In¬ 
clude  (b)  and  (c),  they  cannot  also  Include  (d).  In  fact,  the 
third  competing  collapse  mechanism  is  likely  to  be  a  combination 
of  (c)  and  (f),  because  this  allows  us  to  increase  the  co¬ 
efficient  of  Mq  at  the  expense  of  that  of  Mq  and  again  requires 
a  lairge  dissipation  of  work.  We  try  (g)  »  (c)  -2(f)  (see  Table 
III),  obtaining  the  Inequality 

4Mq  +  2Mq  +  4Mq  -  520.  (16) 


t 
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The  coefficients  in  (15) »  (1^)#  and  (l6)  can  be  combined  with 
the  positive  factors  2,  2/5,  and  8/5  to  yield  the  coefficients 
in  (12).  The  mechanisms  (b),  (c),  and  (g)  thus  are  admissible 
as  competing  collapse  mechanisms .  To  investigate  whether  they 
yield  the  solution  of  our  problem,  we  solve  the  equations  ob¬ 
tained  by  using  the  equality  sign  in  (15),  (l^),  and  (l6). 

f 

TSius, 

Mq  -  58.5,  Mq  -  55.5,  m"  -  15.0  (17) 

(in  ft.t.).  Using  these  values  with  the  appropriate  signs  for 
the  bending  moments  at  the  plastic  hinges  of  the  three  mecha¬ 
nisms,  we  compute  the  remaining  bending  moments  from  the  equa— 
tions  of  equilibrium  obtaining  the  entries  in  row  M  of  Table  III. 
Since  none  of  these  exceeds  the  relevant  fully  plastic  moment 
in  absolute  value,  the  design  (17)  has  the  minimum  structural 
weight . 

We  conclude  this  introduction  to  limit  analysis  and  design 
with  a  brief  history  of  the  field.  A  vague  belief  that  a 
structure  cannot  fall  as  long  as  there  exists  a  distribution 
of  stresses  that  are  in  equilibrium  with  the  applied  loads 

I 

I 

and  do  nowhere  exceed  the  strength  of  the  material  underlies 
the  methods  of  structural  analysis  that  were  in  use  long  before 
the  advent  of  the  theory  of  elasticity.  For  an  unusually 
articulate  application  of  this  principle  the  reader  is  referred 

to  a  paper  by  Fourier  [2]  * .  ; 

* 

The  author  is  indebted:  to  Dr.  George  Dantzlg  of  RAND 
for  this  reference. 
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Wlth  the  development  of  the  mathematical  theory  of  elas¬ 
ticity,  this  principle  became  discredited,  and  it  was  only  in 
the  second  decade  of  this  century  that  Kazinczy  [5]  in  Hungary 
and  Klst  [4]  in  Holland  pointed  out  that  the  ductility  of  struc¬ 
tural  steel  enabled  a  continuous  beam  to  adapt  Itself  to  loads 
far  in  excess  of  those  under  which  the  yield  limit  was  first 
reached  locally.  The  early  theoretical  and  experimental  work 
in  the  new  field  remained  restricted  to  continuous  beams  (see, 
for  instance,  the  book  by  Van  Den  Broek  [5]  >  who  coined  the 
term  "limit  design"  for  what  is  more  accurately  called  limit 
analysis).  Rigid  frames  were  first  investigated  from  this 
point  of  view  by  Baker  [6] ,  whose  full-scale  experiments  con- 
tributed  much  to  the  acceptance  of  the  method  by  structural 
engineers.  The  fundamental  theorems,  reflecting  the  static 
and  kinematic  points  of  viev^  were  first  rigorously  established 
by  Greenberg  and  Frager  [7]  for  beams  and  frames  and  soon 
generalized  to  other  types  of  structures  by  Frager  [8]  .  The 
use  of  a  complete  set  of  linearly  Independent  collapse  mechan¬ 
isms  was  Introduced  by  Neal  and  Symonds  [9]  .  For  comprehen¬ 
sive  '.presentations  of  theory  and  practice  of  limit  analysis 

I 

i 

fran.ithe  point  of  view  of  the  structural  engineer  the  reader 
is  referred  to  the  texts  of  Baker,  Home,  and  Heyman  [lO] , 

Neal  [11] ,  and  Heyman  [l] ,  which  contain  numerous  further 
references . 

The  methods  of  limit  design  are  more  recent  and  less  weir^-. 

I 

developed  than  those  of  l^mlt  analysis.  Heyman  [12]  outlined 
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the  static  approach,  while  Foulkes  [15]  developed  the  kine¬ 
matic  approach  and  recognized  the  relation  to  linear  program¬ 
ming.  Livesley  [14]  proposed  a  method  of  steepest  descent 
towards  the  structure  of  minimum  weight  and  coded  it  for  the 
EDSAC  at  Cambridge  University.  For  a  particular  example, 
Prager  [15]  investigated  the  influence  of  the  linearization 

f 

of  the  expression  for  the  structural  weight. 


p 


P-11 2^ 

7-11-57 

-22- 


i 


P-1125 

7-11-57 

-26- 


BIBLIOGRAPHY 


1.  Heyman,  J.,  Plastic  Design  of  Portal  Frames,  Cambridge 

University  Press,  1957* 

2.  Fourier,  J.  B.  J.,  Oeuvres,  Gauthler-Vl liars,  Paris,  I89O, 

vol.  2,  p.  517. 

5.  Kazlnczy,  G.,  Betonszemele  (Budapest)  2  (1914),  68,  85, 

100. 

4.  Klst,  N.  C.,  Inaugural  Dissertation,  Technical  University 

of  Delft,  1917. 

5.  Van  Den  Broek,  J.  A.,  Theory  of  Limit  Design,  Wiley,  New 

York,  1948. 

6.  Baker,  J.  P.,  Journal  Institution  Civ.  Engrs.  31  (19^9) 

188. 

7.  Greenberg,  H.  J.,  and  W.  Prager,  Proceedings  Am.  Soc.  Civ. 

Engrs.  Separate  No.  59  (1951) »  reprinted  with  discus¬ 
sion  InTrans  act  ions  Am.  Soc.  Civ.  Ehgrs.  117  (1952)  447. 

8.  Prager,  W.,  Sectional  Address,  8th  Intemat.  Congr.  Appl. 

Mech.,  Istanbul,  1952j  Published  In  vol.  2  of  the 
Proceedings  of  the  Congress,  Istanbul,  1956,  p.  65; 
see  also  Id.,  Journal  des  Mathematlques  Pures  et  Appl. 
(1955)  59^7 

t 

9.  Neal,  B.  G.,  and  P.  S.  Symonds,  Journal  Franklin  Institute 

^  (1951)  383. 

10.  Baker,  J.  P.,  M.  R.  Home,  and  J.  Heyinan,  The  Steel  Skeleton 

vol.  2,  Cambridge  University  Press,  1956’! 

11.  Neal,  G.  The  Plastic  Methods  of  Structural  Analysis, 

Chaimian  and  Mall,  London,  19^6. 

12.  Heyman,  J.  Quarterly  of  Appl.  Math.  8  (1950)  373; 

Structural  .Engineer  ^  (1953)  125- 

13.  Poulkes,  J.,  Quarterly  of  Appl.  Math.  10  (1952)  3^7; 

Proceedings  Roy.  Soc.  A  223  (195^) 

14.  Llvesley,  R.  K.,  Quart.  Journal  Mech.  Appl.  Math.  £  (1956) 

257. 

15.  Prager,  W.,  Journal  Eng.  Mech.  Division  Proceedings  Am.  .. 

Soc.  Engrs.  82  (1956)  Paper  1073* 


